Stresses in silos: Comparison between theoretical models and new experiments 
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We present precise and reproducible mean pressure measurements at the bottom of a cylindrical 
granular column. If a constant overload is added, the pressure is linear in overload and nonmonotonic 
in the column height. The results are quantitatively consistent with a local, linear relation between 
stress components, as was recently proposed by some of us. They contradict the simplest classical 
(Janssen) approximation, and may pose a rather severe test of competing models. 

PACS numbers: 46.10.+Z, 83.70.Fn 



The prediction of static stresses in dry, cohesionless 
granular matter has become the focus of renewed atten- 
tion (see 0-^). Surprisingly, there is no consensus on 
what is the basic physics involved. Some argue that the 
behavior is essentially elastic (ultimately justified by the 
slight elastic deformation of individual grains) oth- 
ers that it is dominated by the extremely nonlinear con- 
straint that tensile intergranular forces are absent [^,1). 
Indeed, some of us have argued that the statics 

of granular materials can be described, without consid- 
ering elastic displacements, by assuming a local, history- 
dependent, relation between stress tensor components. 
This gives hyperbolic equations for the stress field, in 
contrast to the elliptic (or elliptic-hyperbolic) equations 
of conventional elastic (or elastoplastic) models. Our 
approach provides a simple continuum model of 'force 
chains' ]9|-pd||; (physical) force chains become (mathe- 
matical) characteristics of the hyperbolic equations. In 
the simplest case, these form a regular array; stresses 
propagate through space via a wave equation ||^,^. Ac- 
cording to the model, the medium is 'fragile' in a precise 
sense j^]: it responds linearly to a specific class of 'com- 
patible' loads; all others cause plastic reorganization. 

This approach accounts well ||] for the pressure 'dip' 
below the apex of a conical sandpilc poured from a point 
source Q. (It also predicts that the dip is absent 
for a pile made of successive horizontal layers, as recently 
confirmed by experiment p^.) However, it has excited 
strong criticism in some quarters 0], and certainly de- 
mands further experimental test |1I|7 For example, such 
models predict that if a small localized overload is placed 
on top of a granular layer, the excess weight at the bot- 
tom is maximal, not directly beneath the weight, but on 
a ring 1^,0. To test this directly is difficult, because of 
strong nonlinearity and (especially) noise effects which 
hinder the interpretation of data 14 1^. 

A more robust and practical situation, is the cylin- 
drical granular column, or bin. Here also noise effects 
come into play; but ways around these (by careful ensem- 
ble averaging of experimental data) have been pioneered 
in pM. Below we report precise measurements (beyond 



those of iQ) of the effective mass Me, supported by the 
bottom plate, as a function of the total mass Mt poured 
into a (small) bin, with and without an added overload. 
With no overload, as expected, Me(Mt) first rises lin- 
early, then saturates at a column height comparable to its 
width; for high bins, most of the mass is 'screened' by fric- 
tional transfer to the walls. A simple hyperbolic model 
(called OSL for 'oriented stress linearity' |^) gives bin 
results close to, but different from, the classical Janssen 
approximation (recalled below) . In contrast to tra- 
ditional methodologies JlSt our new ensemble-averaged 
experiments can distinguish these predictions; we find 
that OSL, which has an extra fitting parameter, is dis- 
cernibly better. Another classical model (ife, see below) 
gives wholly inadequate answers unless unphysical values 
of the wall and bulk friction constants are used. 

There then follow, from the OSL model, two important 
new predictions for the effect of a uniform overload of 
mass Q at the top of the granular column. First, Afg 
should be linear in Q; second, for large Q, Me should be 
nonmonotonic in Mt. We find that, with no further fit- 
ting, our overload experiments quantitatively confirm the 
OSL predictions, strongly supporting the hyperbolic pic- 
ture. At the end of this Letter, we comment on the chal- 
lenge these new results pose to other modelling strategies. 

First we recall our own approach. By stress continuity. 
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where aij is the (symmetric) stress tensor, p is the den- 
sity of the material, and gj is the gravitational accelera- 
tion. In general one needs extra physical assumptions to 
close Eq.0. For an elastic body, one assumes a (single- 
valued) displacement field, and a linear relation between 
stresses and strains (Hooke's law). For poured cohesion- 
less grains, the definition of a macroscopic displacement 
is problematic (see ||l|,Q). Instead we assume that the 
arrangement of granular contacts gives, on continuum 
length scales, a definite relation between components of 
the stress tensor J7|,p0||^ . One such relation, often used in 
the literature, is the ife ('incipient failure everywhere') 
assumption: that the material is everywhere on the verge 
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of Coulombic failure (see e.g. pl|,p|). Then there ex- 
ists a (locally varying) set of axes n ± m such that 
Cnm = <7„„ tan (p where is the Coulomb angle. 

Our modelling strategy instead gives a fundamental 
role to the network of force chains which, if grains are 
undeformable, must carry forces longitudinally One 
interpretation of our equations is that the friction be- 
tween parallel force chains is fully mobilized; a Coulomb- 
like condition, cr„,„ = Onn tani/), then holds (with 'ip < (p 
an 'effective' friction angle) but the orientation m, which 
is directed along the force chains, is now fixed by the con- 
struction history and not (as in ife) by the load P,p7[. 
(This assumes the load is a compatible one.) For simple 
construction histories, like piles and bins, we assume that 
m is the same everywhere, up to an inversion through the 
central symmetry axis; m must then have a fixed angle 
T to the vertical. In cylindrical polars (z, r, 9) with z 
downwards, we recover the OSL model [pi: 
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with 771 = tanr cot(T — ^) and 772 = tanr — cot(T — ?/;). 
Eq.^ closes the problem in two dimensions {d = 2): in- 
serting it into Eq.|l], gives an anisotropic wave equation, 
with one characteristic along m, and another along a 
direction m' at angle t ~ ijj — tt/2. (These can be inter- 
changed without affecting Eq.^; so m' describes a second 
family of force chains |^.) For c? = 3, a further clo- 
sure equation is needed. Our choice here 
but from work on conical piles, we expect insensitivity 
to this choice ppl||. In the bin geometry, the OSL model 
can then be solved exactly {d — 2) or numerically 

(d = 3). Note that ife, like OSL, gives propagative (hy- 
perbolic) equations; but these are nonlinear, unlike our 
wave equation. 

For nonzero 772, the force chain network distinguishes 
between inward and outward radial directions. This does 
not contradict the axial symmetry present j^. But if as 
well the medium is locally symmetric, then 772 = 0; in 
Eq.^, this recovers the model of Ref . . The latter can 
be viewed as a local version of the classical Janssen hy- 
pothesis pT| , ^ . Janssen proposed a constant ratio be- 
tween horizontal and vertical stresses, Urr — Ka^z, but 
neglected altogether their dependence on r. Assuming 
also that friction at the wall is fully mobilized, with a 
friction coefficient tan0„, he found the equation: 



Moo(l-exp[-Aft/Afoo]) 
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with Moo — pD"^ /2K tan (pyj for d — 2, and Moo — 
pirD^ /16K ta,n pw for d = 3; D is the bin diameter. 

We turn now to the experimental procedure, described 
in detail in The bin is a tube of diameter D — 3.8 
cm, filled with beads of glass (density pi, — 2.6 g/cm"^, 
diameter 2 mm) . The bottom comprises a very stiff scale 
plate (2 X 10'' N/m). Initially, the tube is filled with a 
low packing density; this is increased by giving it small 



taps. The bottom plate is then lowered (by a few tens 
of microns) and the effective mass decreases monotoni- 
cally to an asymptotic value; Mg and the mean density 
p are measured. The density is again increased by tap- 
ping, the plate lowered and further measurements taken. 
This entire procedure is done about 30 times - each run 
giving results for the whole range of densities. The mea- 
sured results for show a certain (Mj-dependent) 'error 
bar': not a measurement error of the mass, but arising 
from intrinsic fluctuations in the packing. This proto- 
col is a major advance because (a) an ensemble average 
value for Me is found, improving accuracy; (b) due to 
the downward motion of the base, that wall friction is 
fully mobilized, which might not otherwise be the case 
p3| . The wall friction angle is measured separately as 
pyj — 22° ± 2° p6|, thus eliminating one fit parameter. 
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FIG. 1. Main figure: Experimental and theoretical Me (Mt) 
curves, (a): Statistical dispersion of the measures, (b): Rel- 
ative deviation between experiment and theory, i.e. SM^/A. 

The experimental results, for a packing density p — 
1.53 g/cm'^, are compared in Fig.l with three models: 
IFE (which has no adjustable parameter once the inter- 
nal friction angle p> — 25° ±2° is known); Janssen's equa- 
tion (one adjustable parameter); and the OSL model (two 
adjustable parameters). Each plotted datapoint is itself 
a mean value, with an error bar A shown in inset (a). 
(This is small at small Me but then grows rapidly.) To 
find the best fits, we have minimized the following: 



= A^-i^((5M*/A' 
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where SMI is the difference between the ith experimental 
datapoint and the theoretical Me value. A' the observed 
error bar, and N the number of datapoints. 

For our data, the (active) ife approach, using the mea- 
sured friction values p and p^, is plainly inadequate. Bet- 
ter agreement with ife is found by taking p> and/or (p^, 



2 



as fit parameters. Even then, the fit remains poor (e.g. 
E — A. 4:3 for p = 1.53 g/cm^) ; and the fitted values, 
4' — 4'w — 30° are incompatible with those found by di- 
rect experiment. For given (j)^, ife systematically over- 
predicts the asymptotic stress; so the fitted 4>w exceeds 
the real one. In systems where the wall friction is not 
fully mobilized, the error is harmlessly absorbed by the 
fit. In our system, the fitted value is higher than the fully 
mobilized (p^ measured separately, which is unphysical. 

Unlike the ife model, Janssen's model gives a fair ap- 
proximation {E ^ 2; Table 1) but, as shown in inset 
(b), there is a clear systematic deviation: screening by 
the walls is in turn over- and underestimated for small 
and large Mt values. (Note also that our K parame- 
ters are higher than those usually reported |18 1^: but 
as with IFE, low fitted values might compensate for in- 
completely mobilized of wall friction.) This has led two 
of us to propose elsewhere an empirical model (not 
shown) where an excess contribution from grains at the 
bottom of the pile is added to the Janssen result. As 
shown in Table 1, the best-fit OSL model does as well 
as this empirical model, with an error E ^ 1 [Q: the 
systematic deviations are reduced, in particular in the 
first part of the curve. This can be understood by noting 
that within the OSL model, the grains contained within a 
'light-cone', resting on the bottom plate, cannot interact 
with the walls ; the mass of these grains is completely 
unscreened. 
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p = 1.51 g/cuv^ 
{p/pb = .58) 


E = 5.96 


E = 2.11 
Afoo - 61.9 g 
K==0.65 


E = 0.89 
T]i = 0.55 
m = 1.03 


p = 1.53 g/cm-^ 
{p/pb = .59) 


E = 8.55 


E = 2.28 
Moo - 55.3 g 
K=0.74 


E = 0.94 

771 = 0.64 

772 = 0.97 


p — 1.56 g/cw? 
{p/Pb = .60) 


E = 10.1 


E = 2.28 

= 52.3 g 
K=0.80 


E = 1.02 
7/1 = 0.71 
T]2 = 0.85 


p = 1.59 g/cnT' 
{p/Pb = .61) 


E = VIA 


E = 2.30 
Moo = 48.5 g 
K=0.87 


E = 1.08 

771 = 0.87 

772 0.49 



Table 1. Results of the fits of the experimental data points, 
and the corresponding physical parameters. 

Note the values found for 772. The minimum of E{rj2) 
is not sharp, but positive 772 is always preferred (as for 
other types of grains [|l^). A positive 772 means that 
most of the weight follows the 'inward' characteristic thus 
reducing the screening effect of the walls. Conversely, 
in sandpiles (created from a point source) 772 is negative 
1^,0; this 'outward' transfer of weight is responsible for 
the pressure dip underneath the apex. Positive 772 could 
be caused by slight inward avalanches of material as the 
base is lowered. Its decrease at higher densities might 



indicate a diminished susceptibility to this effect; alter- 
natively the tapping procedure could progressively erase 
a local assymmetry induced by the initial fill. 




OSL model (ti,=0.76 and ii2=0-74) 
OSL model for different overloads 
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FIG. 2. Me(Mt) for various overloads Q. Symbols: exper- 
iment (p = 1.60 g/cm'^); long-dashed lines: OSL predictions. 
Inset: (Me — Moofo)/Q showing linear data collapse. 

We now turn to the key results of this paper, for the 
response to an overload Q placed on top of the granu- 
lar column. (This is a solid piston, just narrower than 
the cylinder.) This is taken into account within the OSL 
model by modifying the boundary conditions to include 
a uniform downward stress at the top surface. Such a 
load is found to be compatible. From the linearity of the 
OSL model (also true of Janssen's model) we then have: 



Me = Moo/c 
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In Janssen's description, /o(a;) = 1 — e~'^ and fQ{x) = 
e~^, so Me is monotonic in the poured mass Mt (and 
constant when Q — M^o)- The result of the OSL model 
are more surprising: /o(a;) and fqix) have different x 
dependences. Hence M,, is not monotonic in Mt; at in- 
termediate Q there is an 'overshoot' (Fig. 2). In addi- 
tion, both functions have a (slight) oscillatory character, 
caused by 'resonances': these are standing- wave modes of 
the wave equation, damped by 'absorption' arising from 
wall friction (see and (l^). In Fig. 2, we show the ex- 
perimental results obtained for various overloads Q. As 
shown in the inset, these results do indeed obey the linear 
relation, Eq.||, to good accuracy. A clear overshoot effect 
is also seen, although any further 'resonant' oscillations 
are small (even theoretically). Note that the OSL pre- 
dictions in Fig. 2 use the same parameters as determined 
previously for Q = 0. Thus OSL, with no further fitting, 
gives a good quantitative account of the data for all Q. 

We have shown that simple hyperbolic models 
encoding the presence of linear force chains 0,S| , can be 
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used to reproduce quantitatively the observed stress re- 
sponse of cohesionless granular media, not only in piles 
1^, but in bins. The same is not true of the traditional 
Janssen analysis. Nor is it true of ife; this does predict 
resonant behaviour (at least in local stresses ||2l[]), but 
our results, even without overload, rule it out entirely as 
a physical model. Any expectation of nonmonotonicity 
in Fig. 2 based on ife would thus have been misplaced. 

What of other continuum modelling strategies? Much 
recent work on bins and silos has studied elastoplastic 
constitutive models (also widespread in soil mechanics), 
often by a finite-element method. There are many such 
models, and a recent comparative study found little con- 
sensus among them p5| ]. But we wonder whether these 
approaches can, with reasonably few fit parameters, re- 
produce the results of Figs.l and 2. For example, the 
observed linearity in Q (seen even for Q/Moa — 1) may 
set a challenge, although one finds numerically that, af- 
ter summing stresses over the base, the (non-linear) ife 
model obeys to a good precision the linear relation (||). 
Linearity is, of course, also recovered if the material is en- 
tirely Hookean. The challenge is then to explain within 
a purely elastic theory the nonmonotonic (if not oscilla- 
tory) curves of Fig. 2. The investigation of these impor- 
tant questions is underway p6| . 

Finally it is important to map out more clearly the do- 
main of validity of the hyperbolic approach (see e.g. pof). 
In particular, our granular columns are tiny: only twenty 
grains or so across. These data clearly do not rule out a 
crossover to more conventional elastic or elastoplastic be- 
havior at larger scales (e.g. where the grains start to de- 
form) |^,|| , although the hyperbolic approach also works 
well in conical piles up to 1 metre wide [ 12| . Careful over- 
load experiments on larger bins could be very valuable, 
as well as local stress measurements, which would reveal 
more clearly the oscillatory nature of the response. 
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